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Abstract 

We obtain multi-soliton solutions with discrete symmetries in the chiral quark soli- 
ton model using the rational map ansatz. The solutions exhibit degenerate bound 
spectra of the quark orbits depending on the background of pion field configurations. 
It is shown that resultant baryon densities inherit the same discrete symmetries as 
the chiral fields. Evaluating the radial component of the baryon density, shell- like 
structure of the valence quark spectra is also observed. 
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1 Introduction 

The chiral quark soliton model (CQSM) was developed in 1980's as a low- 
energy effective theory of QCD. Since it includes the Dirac sea quark con- 
tribution and explicit valence quark degrees of freedom, the model interpo- 
lates between the constituent quark model and the Skyrme model [1,2,3,4,5]. 
The CQSM incorporates the non-perturbative feature of the low-energy QCD, 
spontaneous chiral symmetry breaking. It has been shown that the B = 1 so- 
lution provides correct observables such as a nucleon including mass, electro- 
magnetic value, spin carried by quarks, parton distributions and octet SU(3) 
baryon spectra. For B = 2, the stable axially symmetric soliton solution was 
found in Ref. [6]. The solution exhibits doubly degenerate bound spectrum 
of the quark orbits in the background of axially symmetric chiral fields with 
winding number two. Upon quantization, various dibaryon spectra were ob- 
tained, showing that the quantum numbers of the ground state coincide with 
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those of physical deuteron [7,8]. For B > 2, the Skyrmc model predicts that 
minimum energy solutions have discrete, crystal- like symmetries [9,10,11]. Ac- 
cording to this prediction, we studied the CQSM with B = 3 tetrahedrally 
symmetric chiral fields and obtained triply degenerate spectrum of the quark 
orbits [12]. Its large degeneracy indicates that the tetrahedrally symmetric 
solution may be the lowest-lying configuration. Thus, for B > 3, one can 
also expect that the lowest-lying solutions in the CQSM inherits the discrete 
symmetries predicted in the Skyrme model and have large degeneracy. 

In this paper we shall study sohton solutions with S > 3 in the CQSM using 

the rational map ansatz for the background of pion fields. The procedure of 
the numerical computation performed here is shown in Sec. 3. In Sec. 4, we 
show obtained classical self-consistent soliton solutions with B = 3 — 9. These 
solutions exhibit various degenerate spectra of the quark orbits depending 
on the symmetry of the background chiral fields. Such degeneracy generates 
large shell gaps, which suggests that the solutions are stable local minima. 
Evaluating the radial component of the baryon density, shell-like structure of 
the valence quarks can be observed. The results and discussion are summarized 
in Sec. 5. 



2 The model 



The CQSM incorporates the nonperturbative feature of the low-energy QCD, 
spontaneous chiral symmetry breaking with the vacuum functional defined by 

[1] 



Z = J VTrVijjVilj^ exp i j d*xi}{i^ - MW') 
where the SU(2) matrix 



(1) 



[/TB ^ l±^c/ + ^— ^C/t with U = exp (iT-Tr/U) 

describes chiral fields, t\} is quark fields and M is the constituent quark mass. 
/tt is the pion decay constant and experimentally /^r ~ 93MeV. 

The B = \ soliton solution has been studied in detail at classical and quan- 
tum level in Refs. [1,2,3,4,5]. To obtain solutions with B > 1, we shall em- 
ploy the chiral field configuration with winding number B obtained in the 
Skyrme model as the background of quarks, which can be justified as follows. 
In Eq. (1), performing the functional integral over ip and fields, one obtains 
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the effective action 



-iNcSplniD = -iNclogdetiD, 



(2) 



where iD = i(f) — MU^^ is the Dirac operator. The classical solutions can 
be derived by imposing an extremum condition on the effective action with 
respect to U . For this purpose, let us consider the derivative expansion of the 
action [5,13,14]. Up to quartic terms, we have 



where = Wd^U. Suitably adjusting the coefficients, one can identify the 
first two terms of Eq. (3) with the Skyrme model action. However, the 4th 
order terms tend to distabihze solutions and no stable classical solution can 
be obtained from the above action [13,15]. Nevertheless, because of their sim- 
ilarity, it will be justified to adopt the configurations of the solutions in the 
Skyrme model to chiral fields in the CQSM. 

In the CQSM, the number of valence quark is associated with the baryon 
number such that a soliton with baryon number B consists oi NcX B valence 
quarks. If the quarks are strongly bound inside the soliton, their binding energy 
become large and the valence quarks can not be observed as positive energy 
particles [16,17]. Thus, one gets the picture of the topological soliton model in 
the sense that the baryon number coincides with the winding number of the 
background chiral field when the valence quarks occupy all the levels diving 
into negative energy region. 

Let us rewrite the effective action in Eq. (2) as 



The effective action Seff{U) is ultraviolet divergent and hence must be regu- 
larized. Using the proper-time regularization scheme [18], one can write 




(3) 



^eff ^ -iN,logdet{idt - H{W')), 



(4) 



and introduce the eigenequation 



(5) 
(6) 
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1/A2 

2 ^ J 27r J T 



(7) 



-OO 1/A2 

where T is the Euchdean time separation, Dq and ifo are operators with U — 1. 
At T ^ OO, we have e**^"^*^ ~ g-i-Bseidr xhe total energy is then given by 

E,tatic[U] = E,4U] + Efield[t/] - ^field[t^ = 1], (8) 



where 



vaX 



(9) 



is the valence quark contribution with the valence energy of the i th 
valence quark, and the vacuum part is 



A 



Afield = NcJ2{ ^iE^.)\E^\ + ^ exp 



47r 



E, 
A 



(10) 



with 



1 _j^f^ (E, 



l2'VA 



A is a cutoff parameter evaluated by the condition that the derivative expan- 
sion of Eq. (7) reproduces the pion kinetic term with the correct coefficient. 



f 



2 NM^ 



dr 



An^ J T 

1/A2 



The contribution to the total energy in the absence of the chiral fields {U — 1) 
can be estimated using the eigenstates of the eigenequation. 



Ho = -ia ■ V + (3M . 



(12) 
(13) 
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In the Skyrme model it is known that sohtons with B > 3 have discrete 
crystal-hke symmetries [9]. Therefore, we expect that sohton sohitions of the 
CQSM inherits the same discrete symmetry as skyrmions. However, it is too 
comphcated to perform a numerical computation if one imposes such discrete 
symmetries directly on the chiral fields. Thus Houghton, Manton and Sutcliffe 
proposed remarkable ansatz for the chiral fields, rational map ansatz [11]. 
According to this ansatz, the chiral fields are expressed in a rational map as 

U (r, z) = exp(iF(r)nfl • r) , (14) 
where 

and R{z) is the rational map. The complex coordinate z is given by z — 
tan(^/2)e*'^ via stereographic projection. 

Rational maps are maps from CP{1) to CP{1) (equivalently, from S"^ to S*^) 
classified by winding number. In Ref. [11] Manton et al. showed that B = N 
skyrmions can be well-approximated by rational maps with winding number 
N . The rational map with winding number possesses (2iV + 1) complex 
parameters whose values can be determined by imposing the symmetry of 
the skyrmion. We shall use this ansatz for the background chiral fields in the 
CQSM. Their explicit forms are presented in Appendix A. Since the chiral 
fields in Eq. (14) is parameterized by polar coordinates, one can apply the 
numerical technique developed for 5 = 1 to find solutions with higher B. 
Detailed calculations will be presented in the subsequent section. 

Field equations for the chiral fields can be obtained by demanding that the 
total energy in Eq. (8) be stationary with respect to variation of the profile 
function F(r), 



5F{r) 
which produces 

S{r) sinF(r) = P(r) cosF(r), (15) 
where 

S{r) = N,Y.ME,)+s\^{E,)^^{E^)){^i\^%{\x\-T)\^x) , (16) 
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P(r) = 7VeEK^(^M) + sign{E^)^^iE^)){^i\^^°J'hR ■ t5{\x\ - r)|//) . 

(17) 

The procedure to obtain self-consistent solutions of Eq. (15) is that 1) solve 
the eigenequation in Eq. (6) under an assumed initial profile function Foi^r), 

2) use the resultant eigenfunctions and eigenvalues to calculate S{r) and P{r), 

3) solve Eq. (15) to obtain a new profile function, 4) repeat 1) — 3) until the 
self-consistency is attained. 

The baryon density b{x) is defined by the zeroth component of the baryon 
current [2]; 



b{x) = -^{ipjoip) = by^{x) + b&e\dix), (18) 



where 



-E«ig^(4°^M4°^)4°H^)Vr(^)] • (19) 



To examine the shell structure of the quarks, we evaluate the radial density 
for the ith valence quark p^'^\r) in which the angular degrees of freedom are 
integrated out, via. 



p«(r) ^ jd^j sin Ode bf,,{r, 9, ^) (20) 
with the baryon number 



B 



^^ydrrV'^r). (21) 



3 Numerical technique 



The numerical method that is widely used in this model is based on the expan- 
sion of the Dirac spinor in the appropriate orthogonal basis. However, since 
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the symmetries with B > 3 are discrete and it is hard to find any symmetric 
operator commuting with the hamihonian, we shaU expand the Dirac field in 
the Kahana-Ripka basis as in the case of B = 1 [16]. The Kahana-Ripka basis 
which was originally constructed for diagonalizing the hamiltonian with the 
chiral fields of B — 1 hedgehog ansatz is a plane- wave finite basis. The basis 
is discretized by imposing an appropriate boundary condition on the radial 
wave functions at the radius Tmax chosen to be sufficiently larger than the 
soliton size. The basis is then made finite by including only those states with 
the momentum k as k < k^g^.- The results should be, however, independent 

on Tmax and /Smax- 

The hamiltonian with hedgehog ansatz commutes with the parity and the 
grandspin operator given by 



K = J + t/2 = Z + <r/2 + t/2, 

where j, I are respectively total angular momentum and orbital angular mo- 
mentum. Accordingly, the angular basis can be written as 



I imM) = E Ci:.Jm) I -a,)) I -n) . 



(22) 



J3T3 



mas 



With this angular basis, the normalized eigenstatcs of the free hamiltonian in 
a spherical box with radius rmax can be constructed as follows: 



u 



(1) 

KM 



,.(2) 



U 



(3) 
KM 



U 



(4) 
KM 



"km 



i3K{kr)\{KK + \)KM) 
\AkjK+i{kr)\{K + lK+\)KM) ) 

ijK{kr)\{KK-\)KM) 
-AkjK-i(kr)\(K -IK- \)KM) 

i/\kjK{kr)\{KK+\)KM) 
-jK+i{kr)\{K+lK+l)KM) 

iAk3K{kr)\{KK-\)KM) ^ 
jK-i{kr)\{K-lK-l)KM)^ 

tjK+i{kr)\{K+lK+'^)KM) 
-Ak3K{kr)\{KK+\)KM) 
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,,(2) 



"km 



^KM 



ijK-i{kr)\{K-lK-\)KM)^ 
^k3K{kr)\{KK-\)KM) ) ' 

tAf,jK+i{kr)\{K +1K + \)KM) ^ 
jK{kr)\{KK + \)KM) 

iAkjK-i{kr)\{K-lK-\)KM) 
-3K{kT)\{KK-\)KM) 



with 



2 max 



{jK+i{kr 



-1/2 



and Ak = k/{Ek + M). 



(23) 



(24) 



The momenta are discretized by the boundary condition (^"i^max) = 0. The 
u, V correspond to the ''natural" and "unnatural" components of the basis 
which stand for parity (—1)^ and (—1)-^+^ respectively. 

Let us construct the trial function using the Kahana-Ripka basis to solve the 
eigenequations in Eq. (6), 



K,, 



K 



lim 



E E E["i^W4l(r,^,^) 



K=Q M=-Kj=l 



+PWMy^'M(r,e,^)]. 



(25) 



According to the Raylcigh-Ritz variational method [19], the upper bound of 
the spectrum can be obtained from the secular equation 



= 0, 



where 



(26) 



A] ;r(K'M',KM)^ J d'xu^^^jHu^SM. 

B] ;r{K'M',KM) = J d'xu%,^Hv%, 



(27) 
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Table 1 

A schematic picture of the matrix elements A{K'M' ,KM) for the case oi B = 3 
(and 5 = 5), up to K, K' = 2. S, Pi, Pq and P_i refer to the elements coupled with 
{K,M) = (0,0), (1,1), (1,0) and (1,-1) respectively. Other elements are all 0. 





(0 0) 


(1 1) 


(1 0) 


(1-1) 


(2 2) 


(2 1) 


(2 0) 


(2-1) 


(2-2) 


(0 0) 


s 












S 






(1 1) 




Pi 












Pi 




(1 0) 










Po 








Po 


(1-1) 








P-i 




P-i 








(2 2) 






Po 




Po 








Po 


(2 1) 








P-i 




P-i 








(2 0) 


s 












s 






(2-1) 




Pi 












Pi 




(2-2) 






Po 




Po 








Po 



Table 2 

Matrix elements B^''^{K'M', KM) for B = 3 (and B = 5) up to K, K' = 2. S, Pi, 
Pq and P_i refer to the elements coupled with [K,M) = (0,0), (1,1), (1,0) and 
(1, —1) respectively. Other elements are all 0. 



(0 0) (1 1) (1 0) (1-1) (2 2) (2 1) (2 0) (2-1) (2-2) 



(0 0) 


S 




S 




s 


(1 1) 




Pi 




Pi 




(1 0) 


Po 








Po 


(1-1) 


P-i 








P-i 


(2 2) 


Po 








Po 


(2 1) 


P-I 








P-i 


(2 0) 


S 




s 




s 


(2-1) 




Pi 




Pi 




(2-2) 


Po 








Po 



j\u,u^ C"'" are given by interchanging u and v in Eq. (27). For Kmax oo, 
the spectrum E becomes exact. Eq. (26) can be solved numerically. 

The angular part of the Kahana-Ripka basis consists of the spherical harmon- 
ics, spin and isospin wave functions. Thus, introducing complex basis, one can 
rewrite the chiral fields in Eq. (14) in the isospin space as 
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2R ^ 2R* ^ 1- 

n± and can be expanded by the spherical harmonics 



M^,0)=E E «^MiKM(^,<^), (29) 

aKM = J n{e, ip)Y*{9, (p) sin 0d9d(f) . (30) 

Then we can perform the angular integration analj^ically with help of integral 
formula for the spherical harmonics given by 



(2/i + l)(2/2 + l) 
^ 47r(2/3 + 1) 



l2m2 ■ 



(31) 



Since the chiral fields in Eq. (14) is less symmetric than the B = 1 hedgehog, 
the hamiltonian has no grand spin symmetry. As a result, the states with dif- 
ferent grand spin couple strongly, and level splitting within the K blocks occur. 
In Tables 1, 2 are schematic pictures of the matrix elements A{K'M',KM) 
and B{K'M' , KM) for S = 3. Although the size of the matrix becomes quite 
large, the functional space can be rearranged and reduced in size owing to the 
symmetry of the chiral fields. For 5 = 3, the space is divided with four blocks 
for each parity. 

Numerically, we have to truncate the expansion of the size by cutoff Ajinax and 
-ftTmax- In the tables, each element consists of a 4A;max x 4A;inax size of matrix. 
This block is spanned by 2 x Y1!k^o{2K + l)/4 for 5 = 3 and hence the total 
matrix size is estimated as 



- (i^max(i^max + 2)/2 X Ak^^ . (32) 

The typical parameter values that we employ in our analysis are /Cmax = 24 and 
-f^max = 8 for 5 = 3, giving rather large matrix size 3800 x 3800 approximately 
. In addition, our analysis is based on the self-consistent scheme and requires 
computational time much more. Obviously the situation gets worse for higher 
B. Therefore, in some cases such as B = 6, 17, we oblige to cut the size of 
the matrix, producing a few % uncertainty in results. In the next section we 
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Table 3 

Mass spectra for S = 1 — 9, 17 also for some excited states B = 5*, 9* (in MeV). The 
data for S = 2 are taken from Ref. [7]. The ratio of the mass E'static to x E^^^^^J 
are compared to that of the Skyrme model [11]. 



B 




















-Enfield 


-^static 


-^static 

Ours 


/-^ -^Static 

Skyrme 


1 


173 


















674 


1192 


1.00 


1.00 


2 


173 


173 
















1166 


2204 


0.92 


0.95 


3 


210 


210 


210 














1633 


3522 


0.98 


0.96 


4 


144 


146 


146 


146 












2628 


4378 


0.92 


0.92 


5 


123 


131 


131 


139 


210 










3265 


5467 


0.92 


0.93 


6 


120 


124 


150 


150 


206 


206 








3740 


6603 


0.92 


0.92 


7 


115 


120 


120 


120 


166 


166 


166 






4554 


7478 


0.90 


0.90 


8 


97 


97 


115 


120 


139 


139 


203 


203 




5229 


8565 


0.90 


0.91 


9 


69 


101 


104 


104 


107 


166 


166 


179 


179 


6046 


9573 


0.89 


0.906 


17 


83 


95 


95 


95 


153 


156 


157 


173 


175 












177 


178 


179 


192 


194 


194 


196 


196 




10586 


18650 


0.93 


0.88 


5* 


157 


157 


157 


232 


232 










2874 


5680 


0.95 


1.00 


9* 


99 


105 


105 


121 


142 


142 


210 


210 


210 


5700 


9742 


0.91 


0.91 



will discuss the convergence properties of our numerical computation in more 
detail. 



4 Results and discussions 

4-1 Numerical results 

Let us first show the results of the spectral flow analysis. For convenience we 
shall take 

F{r) = -TT + irr/X for r < X 

— otherwise (33) 

11 
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r [fin] 

Fig. 2. Self-consistent profile functions for B = 3 — 9, 17. 
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Fig. 3. Valence quark spectra for S = 1 — 9, 17. 

as a trial function for the profile function. In Fig. 1 we show the spectral flows 

for S = 3 — 9. As can be seen, the number of B positive energy levels are 
diving into negative energy region and thus we obtain the baryon number B 
soliton solutions. Putting the three quarks so as to be colour blind on each 
valence orbits as well as all on negative energy sea levels, we performed the 
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Fig. 4. Angular averaged baryon number density of ith valence quarks p'^^\r) of 
i3 = 3 — 9, 17, with the occupation number and the eigenvalue (in MeV). 
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Convergence : E 



3900 



. ,E , versus k 

static val max 

320 



Convergence : E . ,E versus K 




10 12 14 16 18 20 

Discrete momentum number 



3640 

3620- 

> 3600 
o 

a 3580 

^; 3560- 

K 

<a 

~ 3540 
E? 3520 







— •— Tntol 






— 0— Valence level 










• 


r"°" = 6[fin] 

max 


ib) 







max 
260 



Max number of the grandspin 



Fig. 5. The convergence properties of the solution with S = 3. 
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Fig. 6. The convergence properties of the solution with B = 5. i^max is fixed to be 
8 and mmax is varied until convergence is attained. 

self-consistent calculations. The profile functions for S = 3 — 9, 17 are plotted 
in Fig. 2. In Table 3 are the results for the valence quark levels as well as the 
vacuum sea contributions. The valence quark spectra show various degenerate 
patterns depending on the background configuration. In Table 3, we also show 
the results of the ratio of the mass ^'static to 5 x E^^^l^, and comparison to 
that of the Skyrme model. The results are qualitatively in agreement. 

The valence quark spectra for various B are shown in Fig. 3. It is interesting 
that the results strongly suggest the existence of shell structure for the valence 
quarks. The spectra show (i) four fold degeneracy of the ground state labeled 
by Q and various degenerate pattern for excited levels labeled by ^i, • " "i 
(ii) a large energy gap between the ground state Q and the first excited level 
Ai. Small dispersions of the spectra are observed in the results. In some cases 
they are caused by the finite size effect of the basis (ex. 5 = 4). Growing the 
size Tmax and increasing the number of the basis, more accurate degeneracy 
will be attained. 
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All our solutions are local minima obtained under the rational map ansatz in 
Eq. (14). Of course there may exist other stable solutions with lower energy 
outside the ansatz. We, however, suspect that the large degeneracy caused 
by the symmetry of the chiral fields would give a strong contribution to the 
minimization of the total energy. 

In Fig. 4 are the results of p^^\r) for B = 3 — 9,17. The behaviour of the 
density near the origin confirms the existence of three shells {Q,Ai,A2)- Q 
behaves hke "5'-wave" and others hke "P-,£)-wave" in a hydrogen-hke atom. 
However most of the densities are nearly on the same surface and very small 
(not zero) near the origin. The plateau in the density observed at the center 
of the nucleus [20] can not be attained in our solutions. Therefore one may 
need to employ the multi-shell ansatz [21] even in the case of light nuclei. 

Our numerical results depend on the basis parameters: radius rmax, the number 
of discretized momenta fcmax, and the maximum value of the grandspin -f^max- 
In Fig. 5 we demonstrate the accuracy of our numerical computation for B = 3 
with Tmax = 6 [fm]. The convergence of the solution and the valence energy 
-^static, Eyax are shown as a function of (a) /cmax with fixed Kmax and (b) K^aax 
with fixed /cmax- This confirms that for K^^x — 8 and k^ax = 20 the solution 
with B — 3 is sufficiently converging. 

As stated above, dispersions appear in the spectra, and in some cases it will 
be eliminated by taking a larger size of the spherical box and increasing the 
number of basis. However, for most of the solutions obtained here, these dis- 
persions can not be ascribed to numerical errors. Let us examine the relation 
of the energy and the soliton size r^ax for B = 5 (sec Fig. 6). For each vahic of 
'"max, we employ sufficiently large number of the basis to attain convergence. 
As can be seen, the dispersion does not disappear even for larger value of rmax- 
We therefore conclude that these dispersions are not due to numerical errors 
nor uncertainty but are the inherent feature of the solutions. 



4-2 Symmetry and the degeneracy of the quarks 

The bunch of valence spectra due to the potential with discrete symmetries 
has been observed in the study of heavier nuclear systems. In Ref. [22], the va- 
lence spectra are highly degenerate because the deformation of the spherically 
symmetric shell produces large shell gaps. Thus the nuclei can be considered 
to be more stable than the spherical one. As discussed in Ref. [11], the group 
theory should predict the level structure of pion fluctuations. However, our 
problem is more complicated due to the presence of quarks. Before discussing 
it in detail, let us show how the shell deformation is related to the degeneracy 
of the spectrum. 
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In general, if an eigenequation given by 



HiP^ = E^iP^ (34) 
is invariant under a symmetric operation R & g, the equation transforms as 



(35) 



Therefore the states {ip^, Rip^} are degenerate in energy with E^. The set of 

dn eigenfunctions {ipi^^jii = l,...,(i/i) belonging to a given eigenvalue i?^ 
will provide the basis for an irreducible representation of the group g of the 
hamiltonian [23]: 



Rtp 



Y.4^Df{R) 



(36) 



The operator R are constructed as follows. If chiral fields have some particular 
point group symmetry i.e., U{x') — G{a)U{x)G{dy {G{d) G SU{2)i, and a 
denotes the matrix of the following spatial rotation aj^), the Dirac equation is 
invariant under the Lorentz transformation 



[x'y — a'^x^ or x' = ax 



(37) 



with 



X 









/ 


1 










:) 


, d = 












V 






\ 












(38) 



accompanying a corresponding iso-rotation 

{i-f^d,, - MW^{x))il:{x) = ^ {i^dl - MU^^{x'))il:'{x') = , (39) 

with 



iP'{x') = {S{d) X G{d))iP{x) , 



(40) 



where S{d) is a 4x4 matrix which is a function of the parameters of the Lorentz 
transformation a, satisfying aj^7^ = S'^j'^S. The operator R corresponding 
to this rotation is thus defined by 
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i;'{x) =def Rip{x) = S{a) X G{a)ip{a-^x) 



(41) 



One can easily check that R commutes with the hamiltonian in Eq. (6). Con- 
structing R for each symmetry of the hamiltonian, one should be able to 
deduce the degeneracy structure of the spectra occurring in the valence level. 
The construction of R in the case of S = 3 is shown in the appendix B, which 
confirms that the valence level has triply degenerate spectra. 

Our numerical results indicate that the winding number strongly couple the 
elements with different K and hence correlated valence spectra occur (see Ta- 
bic 1-2). As can be seen from the operator K of the B = 2, the degeneracy of 
the valence spectra are determined by the shape deformation (symmetry) as 
well as the winding number of the chiral fields [6] . The four-fold degeneracy 
of the lowest states may be ascribed to the chiral symmetry SU (2) l x SU (2)r 
of the hamiltonian. The degenerate structure for B > 3 will be well under- 
stood if symmetric operators of the hamiltonian which consist of the angular 
momentum, spin, isospin and winding number, are explicitly constructed. 



5 Summciry 

In this paper we investigated the multi-soliton solutions in the chiral quark 
soliton model, using the rational map ansatz as a background chiral fields for 
quarks. The chiral fields with multi-winding number have particular discrete 
symmetries and it was shown that the baryon densities inherit the same dis- 
crete symmetries as the chiral fields. For the quark levels we observed various 
degenerate bound spectra depending on the background of chiral field con- 
figurations. Evaluating the radial component of the baryon density, shell-like 
structure of the valence quark spectra was also observed. The group theory 
should predict these level structures resulting from the symmetry of the back- 
ground potential. In fact the degeneracy of the valence spectra are determined 
by the winding number of the chiral fields as well as the shape deformation 
(symmetry) of solitons. The four-fold degeneracy of the lowest states may be 
ascribed to the chiral symmetry SU{2)l x 5'C/(2)r of the hamiltonian. To get 
better understanding of the relation between the quark level structure and the 
winding number or the shape deformation, further analysis will be worth to 
be done in future. 
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A The rational maps 

In this appendix, we present the cxphcit forms of the rational maps which we 
used in our analysis as a background chiral fields (from Refs. [11,24]). 

sfliz^ - 1 

z{z^ - V^t) ' 

+ 2^iz^ + 1 

- 2V3iz^ + 1 ' 
z{z^ + 3Mz^ + 3.07) 
3.07^4 _ 3 94^2 + 1 ' 

z"^ + o.m 

z^(O.WzH + l) ' 

z{z^ + r/Vdz"^ + 7z^ - 7/^5) ' 

z^ - 0.14 
^2(0.14^6 + 1) ' 

z(-3.38- 11.19^^4 + z«) 

l-11.19i^4_ 3 38^8 

z{z^ - 5) 
-5z^ + 1 ' 

5iVSz^ - 9z* + SiVSz^ + 1 - 1.98^2(^6 _ ^^^4 _ ^2 ^ ^^3^ 
z^{-z^ - 3iV^z^ + 9^2 _ 5^73 _ L98z{-i^/3z^ + z^ + i^z^ - 1)) ' 

YIz^^ - 187z^° + 119z^ - 1 
^2(^15 + 119^10 + 187^5 + ■ ^ ■ ) 



B The Lorentz transformation with the chiral fields 

In this appendix, we briefly show the evaluation of the rotation operator R 

for 5 = 3 tetrahedron and also present the results of the transformation law 
for the numerical basis (23). The 5 = 3 tetrahedral soliton is characterized by 
two symmetry operations [11]: x and T^. The Lorentz transformation 
operators S and the operators for the chiral fields {^f, /i} corresponding to the 
symmetry operations [x'Y = oj^a;'^ are given by 



i?3 = 
i?4 = 

Rh — 

Rq — 
i?7 = 

Rs — 

R9 = 
i?* = 

Kg - 
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^2 X ^2 : 



1 ^ 
0-100 
10 
-ly 



Sg = Z7°7^^ , 



U (agx) = gU (x)g^ , g = expf-i-ra] 



1 

10 

1 

Y 1 



5"^ = exp[i|-^((723 + CTgi + (712)] , 



.TT 1 



?7(afex) = hU{x)h^ , /i = exp[-«----^(ri + r2 + rs)] . 



The ^ is defined by the direct product of these rotation operators toghther 
with the inverse spatial rotation for the spinor, that is 



Rg^l^{x) ^ (4 X g)i^{a-'x) . (B.l) 
Rhi^ix) = {Sh X h)'il;{al^x) . (B.2) 



We apply these operators to the Kahana-Ripka basis (f) = {«, v} and finally 
obtain the following transformation law: 



Rgi 


KM - 


= (- 


-i\K-Mi 

J-j (pK-M ■ 






'00 = 


000 • 














% -1 ^ 






Rg 


4>b 




0-1 




06 








^10 0^ 






<t>a = 


1 


(011 


+ 01-l) , 


b '- 


= «01O 



(B.3) 
(B.4) 

(B.5) 



V2 



(011 - 0l_l) 
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( ^ \ 



Rh 



\<Pv ) 



V 






-1 











\ 










1 












-1 



























1 

2 


2 















2 


1 

2 


/ 





(B.6) 



i 1 
0^ = ;^(02i + 02-l) , 0,? = ;^(02i - 02-l) , 

i 1 

0C = (022 - 02-2) , 0« = 02O, 0^; = (022 + 02-2)- 

Thus we confirm that the S = 3 tetrahedron exhibits triply degenerate spec- 
tra. 
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